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I. INTRODUCTION ■ 


Flutter is an aeroelastic dynamic instability where both bending and 
torsional components of a wing result in complex uncontrolled motion. Pas- 
sive flutter suppression can be accomplished by increasing structural rigidity 
which requires additional weight. Recent advances in wing modeling and 
control theory have made active flutter suppression more attractive. Active 
flutter suppression to reduce the aeroelastic response of aircraft structures 
results in substantial weight savings, increased aircraft maneuverability, 
rider comfort and gust alleyiation. 

The state-space equations describing the flutter control problem represent 
the flexible structure, unsteady aerodynamics and actuator dynamics, and are 
usually of high order (~60 states). Modern control synthesis techniques 
applied to active flutter suppression, such as the standard linear quadratic 
Gaussian (LQG) solution, are of the same high order as the plant. Implementation 
requirements demand that the full-order control law be reduced or some how 
approximated. Active flutter control designs using the LQG approach and 
some form of control law model order reduction have been reported in Refs. 

1 to 3. The control laws in Refs. 1 to 3 are all analog. 

In this report, an alternative approach to modern control law synthesis 
is applied to the flutter suppression problem. The objective is to demonstrate 
the application of a digital active flutter control design methodology using 
optimal constrained dynamic compensators. The flutter control design is 

I 

synthesized by minimizing a quadratic performance index defined by a weighted 
sum of mean-square steady state responses and control input as in the LQG 
approach. The difference is that the feedback control law is constrained to 
uge only those measurements, y^^, which are practically available. 



Hk = 


( 1 ) 


2k = + Ik 


( 2 ) 


where Xj^ and are the control, state and measurement noise vectors, 
respectively. The wind gusts which continually excite the wing, the measure- 
ment noise, from the sensors and the state initial condition errors are 
specified as part of the design model. Also Included are the "to be designed" 
compensator dynamics, 


-k+1 ®l^-k ^2^-c,k 


'^c , k — k 


( 3 ) 

( 4 ) 


The compensator dynamics are adjoined to the plant model as a completely 
pontrollable and observable stable system. The compensator controls are 
included in and the perfect measurements of the compensator states (except 
for finite arithmetic) are included in _yj^. 

There are a number of advantages to using the LQG output feedback 
approach combined with dynamic compensation. The order of the compensator 
can be specified by the designer subject to stabilizability constraints. 

For the flutter problem, dynamic compensation appears to be necessary to 
ensure stability at high dynamic pressure. The output feedback approach 
allows the designer to choose the simplest and easiest to implement control 
structure which achieves stability margins and low rms response. At the 
same time, dynamic systems which affect the control design, such as actuator 
dynamics and analog prefilters, can be included in the plant model without 
increasing control structure complexity. 

The quadratic cost function used to represent design objectives is 
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specified in continuous-time along with the flutter model dynamics. The cost 
function and plant model are transformed to a discrete optimization problem 
using the sampled-data regulator approach. The resulting optimal output feed- 
back digital control law is a direct digital design. The sample rate can be 
readily changed and new control gains obtained without having to adjust the 
continuous-time cost function weighting matrices to obtain comparable time- 
domain performance. 

The computation delay that occurs in implementing the digital control 
system can be accommodated in the design by weighting control rate in the 
continuous- time cost function. Control rate weighting causes the optimal 
digital control law to use one-step delayed state information. The time 
frame that transpires after receiving the sampled filtered sensor measure- 
ment can be used to compute the control signal that is to be sent to the 
actuator channel at the beginning of the next time frame. 

The optimal contrained dynamic compensator problem has been studied by 
researchers in Refs. 4 to 7. The necessary conditions for the control law 
to minimize the cost function are presented in Ref. 4 for the continuous 
case and Ref. 6 for the discrete case. Despite the many practical advantages 
of the approach, there are only a few serious studies of the application of 
optimal constrained dynamic compensators for continuous plants and none, that 
the authors are aware of, for discrete systems. There are at least three 
reasons which make the application of optimal output feedback design difficult. 
These difficulties are surmounted for the digital flutter suppression control 
law design presented in this report. 

The most Important reason is the lack of a fast, stable, reliable algor- 
ithm for iteratively solving the output feedback necessary conditions for 
relatively large order plants. Recently, in Ref. 8, |;he authors presented 
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a new algorithm for solving the necessary conditions, by computing a 
sequence of gains in Eq. 1 which converges to an optimal gain. This 
algorithm has been applied without difficulty to the digital flutter con- 
trol law synthesis problem. 

Three studies which determine analog output feedback control laws 
using dynamic compensators and cost functions somewhat similar to the 
one in this report are Refs. 9, 10, 11. Gradient search algorithms are 
used to minimize the cost function. 

A comparison of the algorithm used in this report and a typical 
gradient search algorithm is given in Ref. 8. A comparison of the 
symplified version of the algorithm in this report (i.e. the case when 
a is forced to be 1.0 in Eq. 59) and typical gradient search algorithms 
is given in Ref. 12. Both references demonstrate the benefits offered 
by the algorithm in this report. Gradient search algorithms have the added 
difficulty that the search direction can lead to unstable closed-loop plants, 
a problem not encountered by the algorithm used in this report (for a suffic- 
iently small a) . 

A second reason for optimal output feedback design difficulties, 
as reported in most of the references on optimal dynamic compensators, 
is that simply adjoining compensator dynamics to the plant and not 
weighting compensator controls or states in the cost function is un- 
satisfactory since the optimization process causes the unconstrained com- 
pensator gains to approach large values for continuous systems. To avoid 
this problem. Refs. 4 and 5, for example. Include the compensator gains in 
the cost function using a modified version of compensator control weighting. 

A new approach is pursued for the flutter suppression control law design. The 
compensator states and controls are quadratically weighted in the cost function 
and an error term. 
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^ ^-k ■ ^c^k^ 

is added to the cost function. The error term is an attempt to force the 
compensator states to estimate or observe chosen states of the plant. The 
globally optimal dynamic compensator for the case where the order of the 
compensator is n-5, (number of:l plant states - number of plant observations) 
is a reduced-order observer, (Ref. 7). 

A third reason why optimal constrained dynamic compensator may be 
receiving little attention concerns the lack of any theoretical guarantee 
of good control system properties. Theoretical guarantees of phase and 
gain margins for the Linear Quadratic Regulator (LQR) full-state feedback 
control law are discussed in Ref. 13. The stability margins of the full- 
state feedback LQG approach can recover the stability margins of the LQR 
approach by using an input noise adjustment procedure described in Ref. 14. 

A flutter suppression control law designed using the LQG input noise adjust- 
ment method is presented in Ref. 1. Although a proof exists only for a 
full-order LQG control law. Section IVC numerically shows that improved 
stability margins can also be obtained for the discrete optimal constrained 
dynamic compensator approach. The improved stability margins are obtained 
by increasing the initial condition input covariance. Reference 10, which 
determined optimal constrained dynamic compensators by using a state reduced 
full-order Kalman filter as an initial guess and optimizing the gains in the 
reduced-order system, reached a sirailiar conclusion. 

A. OVERVIEW OF THE REPORT 

The state-space model for flutter control synthesis is discussed in 
Chapter II. The large order full-state evaluation model and the lower 
order design models obtained by residuallzation and balancing are presented. 
A baseline flutter model at q = 8.0 kPas, Mach = 0.9, is used throughout 
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the report. The control optimization problem is constructed in Chapter III. 
The sampled-data regulator approach is used. A continuous cost function and 
plant model are transformed to an equivalent discrete-time cost function 
and model. A single-rate, zero-order hold, digital control law design that 
accommodates computation delays is the desired objective. The steps in the 
algorithm used to solve the necessary conditions to minimize the cost function 
are presented in Chapter III. 

Chapter IV presents the baseline control design parameters and closed- 
loop control law performance. The digital controller has a number of para- 
meters which must be varied to trade off stability properties versus rms 
performance. The effects of varying the sampling rate, the prefilter time 
constant and control structure are investigated in Chapter V. As the 
control law must stabilize the wing for a wide range of flight conditions; 
Chapter VI contains stability properties and rms performance for a constant 
gain design and a gain scheduled design from q = 5.0 to 9.0 kPas. Chapter 
VII presents some recommendations and conclusions. Appendix A investigates 
relationships between compensators in this report and observers. 
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II. STATE-SPACE MODEL FOR FLUTTER CONTROL SYNTHESIS 


The control law s}mthesis method presented in this report is used to 
S3nithesize an active flutter suppression control law for an aeroelastic wind- 
tunnel wing model. The geometry of the sweptback, canilevered wing model, 
along with sensor (accelerometer) and control surface locations, is shown in 
Fig. 1. The half-wing is scaled to flutter within the operational limits of 
the NASA Langley Transonic Wind-Tunnel. The wing has an electro-hydraulic 
servo-actuated trailing edge control surface. The surface hinge line is 
located at 80% of the local streamwise chord. The reaction torques of the 
actuator are constrained by a link to the main structural beam in the control 
pod section. The main structural beam is a single tapered aluminum bar con- 
struction with a cruciform cross section as shown in Fig. 1. The accelerometer 
and control surface locations were provided by NASA as part of the wing model. 

A. AEROELASTIC MODEL OF THE WING 

The mathematical model of wing deformation used in this report is discussed 
in detail in Ref. 15. Previous analyses of the problem have employed general- 
ized coordinates, based on zero airspeed vibration modes or other fixed wing 
deformation shapes, from which generalized aerodynamic forces have been computed 
(Refs. 16-18). The model in Ref. 15 employs physical coordinates of bending and 
torsion of the wing structure directly, and uses constant influence coefficient 
matrices to describe the structural, inertial and aerodynamic forces. 

The model is constrained to 7 node points in 2 degrees-of-freedom (vertical 
deflection and rotation in the flight direction). Two additional states are used 
for aileron deflection and an internal housing state linking the control surface 
to the main structural beam. The aeroelastic effects of the aileron are not ne- 
glected when the model is constructed, while the internal housing state has no 
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aeroelastic effect. The equations of motion are generated in terms of the 
forces and moments affecting the 15 degrees-of-freedom. 

The state space model for the flutter problem can be written in standard 
form as 




16 

16 

15 

T 

*T 

.T 

tI 

X = 

z , 

z > 



— w 

-W 

-qa 


and is a 47x1 state vector. The 16x1 vector, z , contains the 7 vertical 

— w 

deflections, h^, 7 rotations, a^, an internal housing state linking the con- 
trol surface to the main structual beam, 6^, and actuator position, 6^. The 

15x1 vector, x , is the unsteady life and moment vector. The 2x1 distrubance 
— qa 

vector, x^, apd 2x1 Dryden wind model vector, 2ig are states for the external 
disturbance models which affect the unsteady aerodynamics. Equation 7 re- 
presents the accelerometer measurement. The parameter, in Eq. 7 is 
accelerometer measurement disturbances modeled as zero-mean white noise with 


covariance V . The scalar, p, is a zero-mean white noise driving term for 

di 

the Drydep wind model and has unit covariance. 
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The second-order Dryden wind model used to represent the Von Karmen power 


spectrum is 

n- (n- 


Is 

n 






(9) 


a - root-mean-square (ms) gust velocity -0.3048 m/s 
wg 

L - scale of turbulence or characteristic length -30.48 m 
- flight velocity, m/sec 

The actuator dynamics are represented by a third-order transfer function 
given by (Ref. 1): 


X 


214 


s + 214 


( 10 ) 


89450 


X 


s + 179.45s + 89450 


( 11 ) 


The first-order pole in Eq. 10 is modeled as the state x in Eq. 6. The 
second-order polynomial in Eq. 11 is embedded in the wing model state-space 


representation. The parameter A in Eq, 

3 

Eq. 10. The matrices A and B in Eq. 6 

g g 

representation of Eq. 9. 


6 has the value 214.0 as shown in 
are detemined from the state-space 


B. MODEL STATES CAUSED BY THE CONTROL STRUCTURE 


The accelerometer signal used for feedback senses low and high frequency 
wing motion as well as noise. Flutter is predicted at a dynamic pressure of 
5.36 kPa and a frequency near 50 rad/sec (the lowest wing motion frequency 
in the wing model) . Digital control laws can have improved perfomance if a 
properly chosen analog prefilter is used to suppress high frequency distur- 
bances in the sensor output before the sensor output is sampled, (Refs. 19 
and 20) . 
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A first-order prefilter with the transfer function 


y 


f 


a 

s + a ^ 


( 12 ) 


is used in the digital control law design. If a full state feedback linear 
quadratic regulator control system is designed, then all plant states, in- 
cluding prefilter states, would have to be measured for feedback, however, 
the use of output feedback does not require the measurement of all the states 
and does not neglect the prefilter dynamics. The baseline value for the pre- 
filter pole, a, is chosen to be near the flutter frequency, (-50.0 rad/sec). 

The prefilter also serves to suppress stable high frequency wing deformation 
modes which are aliased near the flutter mode frequency after sampling. 

The digital control law uses dynamic compensation to improve closed-loop 
stability and performance. The order of the dynamic compensator is arbitrarily 
chosen. The compensator states are introduced in the model as stable dynamics 
with perfect control and noise free observation, 

z = -3lz + lu (13) 

— — — c 

y = Iz (14) 

•^c — 


A digital control law requires a finite amount of computation time to 
output a control command after measurements are obtained. Neglecting the 
computation delay can have adverse effects, including instability, on closed- 
loop performance. An effective method for accommodating the computation delay 
in the Linear Quadratic Regulator design approach is to weight control rate 
in the quadratic cost function. Ref. 21. The optimization problem is solved 
by including control rate dynamics in the continuous- time wing model. 

u = V (15) 
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Weighting v in the cost function weights control rate. The discrete control 
law that results from using Eq, 15 requires that the control actuator command 
applied at time tj^ be computed during the time interval, using the 

accelerometer measurement obtained at the sampling instant, t^^ 

C. MODEL-ORDER REDUCTION AND THE DESIGN MODEL 

The complete wing modwl has 47 states to describe wing motion, 5 states 
for disturbance and actuator dynamics and 2+0 states for the control law with 
dynamic compensation (1 prefilter state, 1 control rate state, and a compensator 
states) for a total of 54+0 states. The complete wing model is used to eval- 
uate control performance. A block diagram of the wing model and discrete 
control law is shown in Fig. 2. 

The design model is a state-space representation of the wing that is 
smaller in dimension than the evaluation wing model. The design model is used 
in the active flutter control optimization problem to determine the feedback 
gains in the control law. A lower dimension design model reduces the cost of 
control design without significantly compromising closed-loop performance, as 
will be shown. 

The design model is obtained by residualizing the high frequency and very 
stable modes in the 47 state wing model. The residualization procedure begins 
by computing the eigenvalues and eigenvectors of the 47 state wing model. The 
eigenvalues are shown in Table 1 for the baseline (q = 8.0 kPa ) flight con- 
dition. A real transformation matrix, T, is computed by arranging eigenvectors 
columnwise. If an eigenvector is complex, the real part is placed as one 
column and the imaginary part is placed as the next column in T, 


A X. = X.x, 
w-i 1—1 


= j^ReaKx^) 



+ j 


(16) 


ImagCxj^) 


11 


(17 a,b) 



The transformed wing model is block diagonal and results in the model repre- 
sentation of the wing dynamics. The eigenvalues are ordered as shown in 
Table 1, 


T~^A T = 
w 


a^^ b^ 0 0 

■’’1 “i 

0 0 ^2 ^2 

-b2 3 .^ 


• • 


n 


= A 


w 


(18) 


Complex eigenvalues form 2x2 diagonal blocks in A , while real eigenvalues are 

w 

scalars on the diagonal. 

The transformed model wing model, (A^, B^, C^) is partitioned as 

follows 


^1 


X „ 

— w2 
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>t 

N) 1 

ISJ 

1 L 


^2 

+ 

. - 


X + 

a 


wl 


w2 




(19) 




X 1 
— wl 


^2 


+ D X +E x + Y 
w a w -(j a 


( 20 ) 


The x^ states in Eq. 19 are equal to T ^3^. The derivatives of the high 

frequency and very stable states, x are assumed to be zero and are eliminated 

w2 

from the model, yielding. 


X, =A,t X - + B, x +B. X, 

— wl 11 — wl ul a wl 


( 21 ) 


y = C - X T + (D -C 0^0 0 ^® o) X + (E -C «A “^B o) 

■’ wl — wl w w2 22 u2 a w w2 22 w2' ‘a 


( 22 ) 
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The control design model partition is shown in Table 1 and uses all states up 
to and including those which govern the actuator dynamics. The control law 
produces unsatisfactory closed-loop performance if the actuator dynamics are 
residualized or neglected. A second, higher order evaluation model shown in 
Table 1, is used to compute Bode plots and Nyqulst plots using the program 
DIGIKON, Ref. 22. The version of DIGIKON employed has a maximum limit on the 
number of states allowed in the plant model. 

The residualized procedure used for the wing model is also used to 

reduce the 4 state dynamical representation of the disturbances ^nd x^ 

to one state in the design model, x • The models for x , x,, x and x are 

, g g' -d’ a -g 

numerically balanced, using results from Ref. 23, to improve computational 
accuracy for the design and evaluation models. 

The balancing procedure was also investigated as an alternative method 
for reducing the order of wing model. The control designs using the balanced 
reduced-order models were acceptable, but were poor in comparison to the 
control designs obtained using residuallzation. Closed-loop eigenvalue loc- 
ations for the design and evaluation models were more dispersed for the 
balanced reduced-order design models. 

Grouping all the models together, the instantaneous measurements are 
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and the 14 +a order control design model is 



( 24 ) 
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III. ACTIVE FLUTTER SUPPRESSION OPTIMIZATION PROBLEM 


The control design optimization problem for synthesizing an active 
flutter suppression control system is presented in this chapter. The control 
design uses an infinite-time quadratic cost function to represent design object- 
ives. The quadratic weights in the continuous-time cost function are chosen 
by the designer. The cost function and plant model are transformed to a discrete 
optimization model by assuming the control inputs can only change at equally 
spaced sample points in time. The compensator states and plant states are cross- 
weighted in the cost function in an attempt to make the compensator state 
estimate, or "obseirve'', specified wing model modes. The discrete control 
structure is constrained to use feedback of only specified states. The 
necessary conditions the constrained feedback gain must satisfy in order to 
minimize the discrete cost function is detailed. An algorithm for determining 
a gain which satisfies the necessary conditions, resulting in a local minimum, 
is presented. 

A. CONTINUOUS-TIME OPTIMIZATION PROBLEM 


The continuous quadratic cost function begins with the following standard 
quadratic form. 


J 



Ql 0 
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^1 

Xa 
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(25) 


At this state in development, the above cost function is different from previous 
attempts to design dynamic compensators using output feedback because both ^ 
and u are weighted. 
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If Is not weighted In the cost function, the continuous-time optimization 
problem can become singular. The matrix can be zero in the discrete-time 
optimization problem without causing the singularity problem, but the compen- 
sator gains can still become very large. 

A solution to the output feedback optimization problem using Eq. 25 is 
somewhat ambiguous since the compensator states have no intuitive interpretation. 
Optimal designs using the Kalman filter. Ref. 19, or observer. Ref. 7, use 
compensators with states that have clearly defined roles as observers of plant 
states. A discussion of the problem of when is a compensator an observer and 
what does it observe is discussed in Appendix A. The global optimum for the 
optimization problem with compensator states having arbitrarily specified order 
less than n-£ remains unsolved. 

A heuristic approach is taken for the active flutter suppression control 
design. The approach successfully accommodates algorithm convergence difficult- 
ies, provides compensator states with an intuitive interpretation and contrl- 
bures toward improving control law robustness. A new term, represented in Eq. 

5, is added to the cost function as follows, 


J.iE 


/" 
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X ,X X x^ 

-wl a g f 
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X 1 

— wl 
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( 26 ) 
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Q and are chosen to be full rank. The new term has the following proper- 
ties: 

• As Q is increased, the compensator states are interpreted as more 

closely following (observing) can be arbitrarily chosen 

subject to the rank condition. 

• If Q is nonzero, the algorithm to be discussed in Section III-D con- 
verges to a control law where the compensator and plant are coupled, 
even if the starting stabilizing gain is uncoupled. When Q is zero, 
an uncoupled starting stabilizing gain (K^, K^, and are zero 
matrices in Eq. 31) causes the algorithm to converge to an uncoupled 
gain which satisfies the necessary conditions. 

• The solution with and without Q is non-unic(ue when compensator states 
are included iq the plant model. Different starting gains converge 
to different local minimums. The different compensator designs at 
local minimums are not necessarily related by a transformation matrix 
(Ref. 7). 

• A simple root locus using the wing model at 8.0 kPa ( the design flight 
condition) demonstrated that the wing model at 8.0 kPa cannot be 
stabilized using only the accelerometer measurement. The plant and 
compensator must be coupled to minimize the cost function for the 
active flutter suppression problem. 

The new term in the cost function is similiar to techniques used in 
explicit model following approaches. Ref. 24 and 25. The difference 
is that the compensator dynamics (i.e. the model dynamics in explicit 
model following parlance) are allowed to be changed by the optimiza- 
tion process. 

The continuous problem is transformed to a discrete problem using the sampled- 
data regulator as discussed in the next section. 
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B. DISCRETE OPTIMIZATION PROBLEM 


The controls v and u are assumed to be constant ov^r the fixed 
sampling interval, At. The piecewise constant controls allow the design 
model shown in Eq. 23, and the cost function shown in Eq. 26, to be trans- 
formed to an equivalent sampled-data regulator problem, (Ref. 26), 
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(29) 


J = lim J 
N^ . 


n 


(30) 


The state Eq. 27, as shown in Eq. 23, is equal to the state vector 

l^x^^ x^ Xg upper right partition of F in Eq. 27 is assumed to 

be zero. If this assumption is not used, the implied continuous-time implemen 
tation for u(t) would be a triangular data hold. The control, u(t), is con- 
stant over the sampling instant if F is modified as shown. A white noise term 
with covariance, V, is added to the discrete measurements and is treated 
as a design parameter. 

The class of control laws considered are restricted to be of the form 


18 



C K. K 

y $ u 

L K K 

f c V 


Defining the variables 


6^1 + 


= I + AtK 


K = AtK 
c c 


The Implementable form for the contrpl law with the compensator is 


^+1 ■ '’c >'f.k ''u 

Sk+1 * *„ ''f ^t.k *'c ^ 

C. NECESSARY CONDITIONS FOR OPTIMALITY 

y ' — i ■ - .■ii . ^ ■■ ..I. . I I I • 

One final problem must be resolved before the necessary conditions for 
optimality for the output feedback probelm can be determined. The cost in 
Eq. 29 minimizes the average long term stochastic performance of the plant. 
State and control initial condition errors are averaged to zero as N Increases 
and do not affect the control design. In contrast to stochastic output feed- 
back, the LQG problem and the deterministic optimal output feedback approach 
in Ref. 4, are primarily concerned with driving initial conditions errors to 
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zero; another desirable control objective. A final modification to the 
quadratic cost function is constructed where minimum stochastic performance 
and fast response (and as is shown later, robustness) can be traded off against 
each other. 


, k«l ,k-i-l 

-k ■ ^CL ^ iSo ^CL ^^w % 

i^k = ^k = - ^ ^ ^tk - ^k 

J(K) = -| (J. + J^) 
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(40 a) 
(40 b) 

(41) 

(42) 

(43) 


The plant response is separated into the transient component, arid the 

stochastic component, 2Egj^> shown in Eq. 40 a. Tbe feedback control law 

is similarly partitioned into the transient, stochastic, 

components. is the transient cost with noise sources set to zero while 

J is the average stochastic cost. The objective is to determine the output 
s 

feedback gain K in Eq. 40 b which minimizes the average cost shown in Eq. 41. 
The necessary conditions for minimizing J(K) are straightforward if the 
quadratic weights in and are equal. The tradeoff between and 
is accomplished by changing X^, the covariance of the plant state initial 
condition. 

For the plant dynamics shown in Eq. 23, the feedback gain constraint shown 
in Eq. 31, and the following conditions. 


. 0.0 

Efrw^wTr '^1 = w 6 

|_ w-k -j wj kj 

(44 a, b) 

E . 0.0 

E [4 x[] - 

(45 a, b) 
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(46) 
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The cost in Eq. 41 can be rewritten as follows, (Ref. 8), 


(47) 


J(K) = Js tr 



+ % tr |K’’^(r^ 


p r + 



(48) 


The matrix, P, satisfies the Riccati-like equation, 

P = + Q - ikc - 

UIj v/1j 


(49) 


is the stable closed-loop plant matrix, 

$CL = ^ - TKC (50) 

The matrices $, F and C are defined in Eq. 27. 

The necessary conditions for J(K) to have a minimum are derived in Ref. 
8. There must exist a gain, K, so that is stable. The gain, K, must 
satisfy Eq. 49 and the following 


® ='0> + rKVKV 

j^r"^pr + rJ K fcsc"^ + vj = 


SC 


T 


(51) 

(52) 


A few comments are 

• If W and V are zero matrices, the necessary conditions are the 
solution to the deterministic output feedback problem, J^, with 
zero. 

• If Xp is a zero matrix, the necessary conditions are the solution 

to the stochastic output feedback problem, J , with J zero, and S 

' s t 

r T 

is the covariance matrix for E Xj^ 
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• The covariance, Xq, performs the same function as adding pseudonoise 
to the plant to improve closed-loop performance. 

An algorithm to compute the optimal K which satisfies Eqs. 49, 51 and 52 is 
given in Ref. 8. An optimal gain exists provided an initial stabilizing gain 
exists, r and C have full rank, m<n, £<n and 


W > 0 


Q > 0 


(53) 

(54) 


The above conditions guarantee convergence. The algorithm may still converge 
if the conditions are not satisfied as in the case in this report. 

D. A PRACTICAL CONVERGENT ALGORITHM FOR DISCRETE OPTIMAL OUTPUT FEEDBACK 
Step 1; Choose so that $ - TK^C is stable, aQ£(0,l], 

z > 1, j an integer S 1 and set i = 0. 

Step 2: Solve Eq. 51 for S using K. The Bartels-Stewart algorithm 

in Ref. 27 is recommended. 

Step 3; Solve Eq. 49 for P using 
Step 4; Invert the symmetric matrices: 


p = r pr + R 


(55) 


S = CSC + V 


(56) 


using Cholesky decomposition. If the symmetric matrices P and 

/s 

S are not positive definite, go to Step 8. 

Step 5t Compute d(Kj^) 


■S)EW + 

- hm - h 


T'' -1 
SC^S 


(57) 

(58) 
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Step 6: Compute 


■'i+i - ■'i 


( 59 ) 


Step 7; Evaluate the cost function, J^., using Eq. ^8. If i = 0, set i to 1, 

and go to Step 2. If is negative, go to Step 8. 

If is negative, go to Step 9, otherwise, go to Step 8. 

Step 8: Decrease using z; ~ ot^/z 

Go back to a previous stabilizing gain 


Ki = K^_^, d(K^ = d(K^_j) 


(60 a,b) 


Compute 


■'l+l = ■'l “l"®!) 


Set “ *^1’ ^ 2. 


9 J 

step 9; Compute 


||(K^) * PK^S - [r'^p<j> + SC^ 


If 


(61) 


(62) 


9 J 

^(K.) II and (J, - J. ,)/J, are less than some convergence 
oK. 1 II i X~1 1 


criterion stop, otherwise set ~ ^ ~ 

Step 2. 

Reference 8 derived the property that there exists an 0 < a ^ 1 such that the 
algorithm is stable. The algorithm uses a number of checks to determine if 
the current value of a is too large and reduces a accordingly. Numerical 
experience indicates that the algorithm converges faster for values of a 
closer to 1.0. 
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IV. DIGITAL CONTROL DESIGN 


This chapter presents the design values for the baseline control system. 
The digital control design is determined by choosing alements in Q, R, W, 
and V, finding an initial stabilizing gain and computing the locally optimum, 
K, using the design wing model. The control design is evaluated using the 
evaluation wing model to obtain rms response for states and controls. Nyquist 
plots are determined using a lower dimension evaluation model. Adjustments 
.are made to design elements until a desired performance tradeoff is obtained. 
A. PROCEDURE FOR COMPUTING AN INITIAL STABILIZING GAIN 

The dynamic compensator makes it difficult to apply output feedback 
stabilization procedures to the flutter control problem. A straightforward 
approach is presented that successfully determines an initial stable gain for 
the flutter problem. 

Taking advantage of the block diagonal form of the flutter model, the 
unstable flutter mode in the model is forced stable by adjusting the 2x2 block 
parameters. in Eq. 31 is easily chosen to stabilize and 3 in Eq. 13 is 
set to 10. With all gains but zero, is stable and a new preliminary 
gain can be determined with the algorithm and the modified wing model. 

The process is repeated using the new preliminary stabilizing gain and 
perturbing the flutter mode in the direction of the true value. If the con- 
trol design is relatively insensistive to the artifically Induced plant para- 
meter variations, a valid stabilizing gain can be iteratively computed in a 
few steps. Three iterations were used for the wing model. 

An alternative approach is to start the design at a low dynamic pressure 
flight condition where the flutter mode is stable. Increasing q Iteratively 
and using the previous gain from the algorithm for each increase in q also 
results in a stabilizing gain at the desired flight condition. 
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B. DESIGN VARIABLES AND OBJECTIVES 

The optimal output feedback design problem for flutter suppression has a 
number of design variables which must be chosen by the designer. The wing 
model is similar to (but not the same as) the wing model in Refs. 1 and 10. 

The control performances in Refs. 1 and 10 form reasonable objectives for 
specifying the control design variables and are shox«i in Table 2. The dif- 
ferences between the approaches are 

• Digital design in this report verses analog designs in Refs. 1 
and 10 . 

• 4 compensator + 1 prefilter + 1 control state in this report 
versus 4 compensator + 1 prefilter state in Ref. 1 versus 4 
compensator states in Ref. 10, 

• This report. Ref. 1 and Ref. 10 each use different plant model 
orders for the design wj.ng models. 

The Dryden wind model shown in Eq. 9 disturbs the dynamics for both references. 

The parameter design variables are the prefilter time constant, the 
sampling time At, the accelerometer measurement noise covariance, the cost 
function weights, Q^, Q^, Q, R, R^, R^ and the state initial condition matrix, 
X^. Structural design variables are the order of the compensator, the compen- 
sator observation matrix, H^, and the timing of the y^ observation. 

C. DESIGN VALUES 

Using Refs. 1 and 10 as guidelines, the order of the compensator, a, is 
chosen to be four. The matrix, H^, is chosen so that the first two compensator 
states "observe" the unstable flutter mode in Table 1 and the third and fourth 
compensator states "observe" the second mode in Table 1. The (4 x 10) matrix, 

H^, is chosen with the plant states expressed in reduced-order model coordinates 
as shown in Eq. 26. Elements in are simply zeroes and ones. 
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As a general rule, if poles and zeroes are "close", the pole is usually 
"close" to being uncontrollable or unobservable. Modes 3 and 4, in Table 1, 
represent stable modes that are onJ.y moderately affected by the control design 
as shown in Table 3. Increasing the order of the compensator to account for 
modes 3 and 4 is not required. 

The sampling rate is chosen to be 200 samples/sec. Plant modes with 

natural frequencies above 628.3 rad/sec are aliased by the sampler. The 

aliased frequencies for these modes are shown in the second column of Table 

3, The mapped poles in Table 3 are determined by first computing $ for 

CL* 

the full order closed-loop evaluation model, pomputing the discrete closed- 
loop poles, then mapping the discrete poles back into the s-domain using 
the natural log. 

Discrete pole X = a + jb (63) 

Discrete pole X =Va^ + b^ (64) 

Mapped Discrete pole ^ + ^ ) + j (65) 

Covariance values and quadratic weights chosen for the design are 
shown in Table 4. The covariance, X^, is the input control position ini- 
tial condition covariance. The covariance, X^, is the eleventh diagonal 

element in the X covariance matrix shown in Eq. 47. The covariance X 
o u 

is defined as X^ = E | | . Table 5 shows that increasing X^ Increases 

gain margins and phase margins but also increases control surface activity. 

A compromise must be reached between the conflicting objectives of good 
stability margins and low control surface activity. If the control surface 
pctlvity is high, moderate gusts can cause the surface to saturate quickly, 
the design may become sensitive to unmodeled accelerometer noise disturbances 
and the control system may become susceptible to limit cycling. The design 
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for = 0.10 in Table 5 is selected as the compromise value. The gain 
margins, phase margins and covariance activity in Table 5 compare favorably 
to the analog designs in Table 2. 

Table 5 shows that the input noise adjustment procedure is capable of 
improving stability margins for the output feedback compensator design 
synthesis technique used in this report for the wing model. Improving stab- 
ility margins increases control surface activity. Adjusting the input noise 
affects both the compensator and the output feedback gains. The implication 
is that artificial input noise variation may benefically affect stability 
margins for any design using the output feedback design approach discussed 
in Chapter III (with or without compensator states) . 

The Nyquist diagrams for the three control laws in Table 5 are presented 
in Fig. 3 and show the progressive Increase in phase and gain margins. A Bode 
plot of the reduced-order plant plus the = 1.0 control law design is pre- 
sented in Fig. 4. The rapid changes at 415 rad/sec and 163 rad/sec are caused 
by the zeros near these frequencies as shown in Table 1. The bandwidth in 
Table 5 is the highest frequency in the Bode plot at which the gain is greater 
than -6B. 

Table 6 shows the variance values for deflection and rotation of the 7 
mass bodies which comprise the wing model for a 0.305 m/s (1 ft.) gust input. 
The variances increase from the fuselage to the wing tip but have acceptable 
values. 

D. NUMERICAL EXPERIENCE 

Numerical experience with the algorithm and the flutter control problem 
with dynamic compensation verified that the baseline control design perfor- 
mance given in Table 5 is probably a local minimum. Different starting gains 
converge to 4ifferent designs with different rms performance. The design in 
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Table 5 could probably be improved by starting the algorithm at different 
stable gains and comparing performance. Theoretical and numerical procedures 
which guarantee uniqueness is an area of investigation recommended in Chapter 
VII. 

Figure 5 shows a typical convergence pattern of the output feedback 
algorithm. The largest Improvements in the cost function occur within 10 
iterations. Values for the algorithm's a parameter, defined in Section 
III-D, vary between 0.2 and 0.05 depending on the weighting and covariance 
values. 

In the next section, design parameters are varied and the effect on 
performance is tabulated. Whenever possible, the starting gain is fixed at 
the same value for each variation. 
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V. EFFECT OF VARYING CONTROL PARAMETERS AND STRUCTURE 

In this chapter, the effects of lowering the sample rate, varying the 
prefilter time constant and altering the control structure from a "one-step 
prediction" to an "update" ' implementation are investigated. The starting 
gain column in the Tables indicates if the starting stabilizing gain is 
the same as the starting gain in the baseline design. If the gains are not 
the same, the starting stabilizing gain is the locally optimum gain obtained 
in the previous variation being investigated. 

A, SAMPLE RATE VARIATION 

Slower sample rates lower the requirements for computer specifications, 
computer code efficiency and A/Dj D/A converters. The effect of progressively 
lowering the sample rate from 200 samples/sec to 100 samples/sec is shown in 
Table 7. The complexities of making meaningfvil comparisons, as well as the 
optimization's propensity for advantageous reconfiguration is evident in 
the Table. The control rms response improves with lower sample rates for 
fixed X^, but gain and phase margins degrade, particularly at the higher 
frequency crossover points. The covariance is increased to 0.25 for the 

160 sample/sec design in an attempt to match 6 ^ ^ response with the 200 sample/ 
sec design. The degradation in high frequency gain and phase margin at the 
lower sample rate is evident. Sample rates below 140 samples/sec appear to 
have questipnable performance. 

B. PREFILTER POLE VARIATIONS 

Analog flutter control laws in Ref. 28 used first order analog prefilters 
in all the designs reported. Values for the prefilter pole ranged from -5.0 
to -20.0. The control laws in Ref. 28use two sensors for feedback. Torsional 
wing motion can be made better observable with two sensors. 
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Lower prefilter poles for the one accelerometer sensor used In this 
report degrade rms response for the control surface, but Improve gain and 
phase margins (except for the high frequency gain margin) as shown in 
Table 8. The best compromise appears to be to reduce the magnitude of the 
prefilter pole until response is marginally acceptable. Although not 

studied in this application, a wash-out filter (s/s+b) should be incorporated 
with the prefilter to suppress low frequency motion of the wing caused by 
aircraft maneuvers. 


C. ALTERNATE CONTROL STRUCTURE 

The control law shown in Fig. 2 is purposely designed so that _Uj^ uses 
the Zf feedback. The control Uj^, by definition, is to 

be written by the onboard computer to the actuator output port at the time 
Instant tj^. A one sample time period. At, occurs in real time in the onboard 
computer between receiving Zf from the measurement input port and writing 

Uj^ to the actuator output port. The At time period can be used to compute 
u^ using Eqs. 38 and 39. 

Consider changing Eq. 28 so that k+1^® assumed to be the measurement 
at t^. 
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Equation 66 is modified using the plant dynamics 
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Note that the "measurement noise" and the process noise become cross- 
correlated. The implementable control law which minimizes the cost function 
using Eq. 66 as the measurement vector is 


z, = $ z - + r y, , + r u, , 

— k c — k-1 c f,k u k-1 


'^k = \ vi + yf,k + ^c ^k-1 


( 68 ) 

(69) 


A block diagram is shown in Fig. 6. All multiplications and additions in 
Eq. 69 can be performed before tj^, except the single multiplication 
and the addition to form Uj^. The computation delay is very small. The 
compensator in the control* law in Eq. 68 resembles the "update" Kalman 
filter form. 

The covariance, X^, for the alternate control structure is reduced to 
0.01 to make the control surface variance slmiliar to the baseline design 

control surface variance as shown in Table 9. The reduction in X de- 

u 

grades the stability properties of the alternate control structure design as 
shown by the Bode plot in Fig. 7. The performance of the alternate control 
structure could probably be significantly Improved by adjusting the quad- 
ratic weights. The quadratic weights are not changed in any of the tests 
from those used to optimize the baseline design. The low frequency gain in 
Fig. 7 is unacceptably high, but could be improved using the (s/s+b) wash- 
out filter. 
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VI. EFFECT OF VARYING FLIGHT CONDITION 


The active flutter suppression control system should stabilize the wing 


over a large range of flight conditions. Two approaches are investigated 


for flutter suppression from <j = S.OkPa to 9.5kPa at constant Mach number. 

The first approach uses the constant gain baseline design at q = S.OkPas 
and investigates closed-loop performance with changing flight condition. 

In the second approach, optimal output feedback designs are obtained at four 
flight conditions. A new gain scheduling procedure is used to schedule control 
gains as a function of q by minimizing the weighted difference of optimal 


and gain scheduled closed-loop eigenvalues. 


A. CONSTANT GAIN DESIGN 


The baseline constant gain design at q = S.OkPas, for the design values shown 


in Table 4, has the following values for the block diagram shown in Fig. 2, 
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0.746 0.0994 
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. 0.00942 0.0964 
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(73) 

(74) 
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The open-loop wing model eigenvalues for the first seven modes is shown 
in Fig. 8 . The wing model is stable for q = 5kPa. The closed-loop wing 
model eigenvalue root locuses for the constant gain design are shown in 
Fig. 9, along with the root locus of one of the clqsed-loop compensator 
poles. The constant gain design stabilizes the wing up to q = 9.5kPa 
where the closed-loop compensator pole shown goes unstable. Rms response 
for the control surface, and gain and phase margins are shown in Table 10. 

The rms response remains at a high level as q decreases but the gain and 
phase margins Improve (except for the high frequency phase margin at q = 
S.OkPas) . 

B. A NEW GAIN SCHEDULING PROCEDURE 

A practical method for adapting a control law to changing flight 
conditions is gain scheduling. Control designs are obtained at a number of 
flight conditions which span the operating range of the plant. The gains 
in the control law are scheduled using regression analysis. Flight condition 
parameters that can be measured or estimated in flight are treated as independ- 
ent variables in the regression. 

The previous section demonstrated that a constant gain design stabilized 
the wing over a wide range of dynamic pressure variations. A gain scheduling 
function of the form 

V ■ s S 

should be able to Improve both performance and the stability region. A 
standard regression analysis cost function for computing Gj^ and G 2 is 

J 3 - j! (77) 

The gains, K^, are the optimal output feedback gains obtained at flight 
conditions. is a positive definite symmetric weighting matrix, and tr 
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signifies the trace of a matrix. 

The gains and were computed using Eq. 77 and the identity 
matrix. The closed-loop plant using was unstable within the region of 
the gain schedule. Adjusting the diagonal elements of to stabilize the 
gain-scheduled closed-loop eigenvalues proved to be difficult. An alternative 
method is developed for choosing Q^. The new cost function attempts to match 
closed- loop eigenvalues and eigenvectors instead of matching gain variations. 

Consider two feedback gains, K and K which have the same dimension but 

gs 

are not equal. The closed-loop eigenvalues, Aj , and eigenvectors, x ^ , for K 
satisfy 

($ + rKC) Xj = Aj Xj j = 1 . . . n (78) 

Following concepts in Ref. 29, assume that K and K are related by 


KC X. = K C X. 
-J gs -J 


(79) 


for the eigenvector x ^ . Substituting Eq. 79 into Eq« 78, it follows that K 
and Kgg have the same closed-loop eigenvalue and eigenvector for the plant 

D. 

The cost function in Eq. 77 is altered using the closed-loop eigenvectors 
for the known gains, K^, 


J - .Et tr i(K - K.)(C. X, W. X* cT) (K - K \ 

3 x=l I gs i' ' 1 i i 1 i' ' gs if 


(80) 


The matrix, X^, in Eq. 80 is the closed-loop eigenvectors for arranged 
columnwise. is a diagonal weighting matrix with nonzero positive elements 

along the diagonal. The notation * means take the transpose of the complex 
conjugate of the matrix. The matrix is real if diagonal weights in 
weight an eigenvector and its complex conjugate equally. The cost function 
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minimizes the error between K C x. and C x. instead of between K and 

gs -j i -j gs 

K^. Minimizing the cost function yields 


[<=1 s] 

ll 5^1 

^ - 
’A 


’ N. _ 
lii 


1^1 ^I'^i 

_ 

Nf 

lii 

m 


Nf 

isi 

I. ** 


Inverting the square matrix in Eq. 81 determines the least square solution 
for and G^. If a gain scheduled closed-loop eigenvalue is unsatisfactory, 
the diagonal weight in can be increased in an attempt to force the eigen- 
value to its design value. This gain scheduling procedure is particularly 
applicable to the flutter problem since the primary concern is the stability 
of one unstable mode. 

C. GAIN SCHEDULE DESIGN 

The performances for optimal output feedback designs at q = 8.0, 7.0, 

6.0 and 5.0 kPa are shown in Table 11. The second row in Table 11 shows' 
that rms response can be reduce:d at lower q while maintaining phase and gain 
margins comparable to q = S.OkPat The baseline design and the last three 
designs in Table 11 are used to compute G^ and G 2 in Eq. 81. 

The scheduled gain performance and closed-loop eigenvalue root-locus are 
shown in Table 12 and Fig. 10, respectively. Phase and gain margins are high 
in Table 12 but the resulting rms response is also high. The scheduled gain 
closed-loop plant goes unstable immediately outside the region of flight con- 
ditions used in the gain schedule. The instability is caused by one of the 
closed-loop compensator poles. New gain scheduling procedures is an area 
recommended for continued theoretical Investigations. 
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1711. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 

The conclusions of this study are 

• An efficient stable algorithm for solving the stochastic infinite- 
time discrete optimal output feedback problem can be applied to active 
flutter suppression control design. 

• Practical discrete low-order dynamic compensators can be designed 
using optimal output feedback. 

• A digital control law which accommodates computation delay can 
stabilize the wing with reasonable nils performance and adequate, 
adjustable, galnrand phase margins. 

• The sampling rate for the control law should be at least 140 samples/ 
sec. 

• Lowering the analog prefilter pole increases rms response but also 
improves one of the gain margins and both phase margin?. 

• A new gain scheduling procedure is developed to yield a stable linear 
gain schedule as a function of dynamic pressure over the design flight 
conditions. Despite these developments, further Improvements and 
modifications in gain scheduling are needed since the gain schedule 
did not perform as well as a constant gain design Inside and outside 
the design flight conditions. 

• Designing the flutter control system using a reduced order model, and 
verifying performance of the control design with a higher order more 
exact model, reduced computer cost without significantly compromising 


control performance. 



B. RECOMMENDATIONS 


• Further develop the baseline control designs for both structures in 
Figs. 6 and 2. The goal is to lower control rms response to 4.5 
deg and improve the low frequency phase margin. Avenues for further 
development for the baseline design are 

A Vary the quadratic weights. 

A Increase the prefilter pole to -60.0. 

A Lower the sample rate to 150 samples/sec. 

A Start the algorithm at many different starting gains to 
determine a lower value for the cost function. 

A Vary the initial condition covariances. 

A Try using a dynamic compensator of order two as in Ref. 11. 

A Introduce a second accelerometer measurement as in Refs. 1 
and 28. 

A Move the compensator pole that goes unstable in Fig. 9 
further into the left half complex plane. 

• Test the active flutter suppression digital control in a wind tunnel. 

• Investigate theoretical and numerical developments which may cause 
the optimal dynamic compensator to be unique (up to a similarity 
transformation) . One avenue is to restrict compensator gains so that 
the compensator is an observer as discussed in Appendix A. References 
9 to 11 did not treat all the control law gains as free parameters. 

• Further improve the gain schedule procedure in Section VII-B. Modify 
the approach so that a select number of gains are scheduled using 

Eq. 76 and others are chosen to remain constant. Introduce eigenvalue 
sensitivity into the cost function. 
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A potentially better approach to gain scheduling is to bring the 
selection of Gj^ and in Eq. 76 directly into the control design 
optimization process. ICS has recently developed a multi-model 
approach to stochastic optimal output feedback which could perform 


the optimization. 



APPENDIX A 


COMPENSATORS AND OBSERVER THEORY 


The globally optimal output feedback constrained dynamic compensator 
solution Is currently unresolved If the order of the compensator Is pre- 
speclfled to be less than n-£. If the order Is Increased to n-Z, where Z 
Is the number of observations, then Ref. 7 shows the globally optimum 
dynamic compensator Is an observer which estimates the n-£ state functional 
not measured with the Z measurements. If the order Is Increased to n, 
the optimal solution uses a Kalman filter for the compensator. These 
results suggest the optimal compensator which minimizes stochastic per- 
formance may be some type of observer of a linear function of plant states 
not measured. The cross weighting used In Eq. 26 Is a straightforward at- 
tempt to cause this behavior to occur. This appendix addresses three quest- 
ions about compensators that resulted from the cross-weighting approach: 

• When Is a compensator an observer of a linear function of plant 
states? 

• If the compensator Is an observer, what does It observe? 

• When are two Internal compensator representations related by a 
state transformation and what Is the transformation matrix? 

These questions are resolved using the theory of feedforward control 
developed In Ref. 30. Given a plant 

(82) 
(83) 
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to follow 


the plant trajectory, and control, ^ which causes 
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The feedforward matrices in Eq. 86, if they exist, satisfy the algebraic 
equation 
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Si 3 ... 


Sii$ S^^F S^2 

1 

0 

1 

> 4-1 


-"21 

S 22 

S23 ...^ 


H D 0 . . ._ 


( 87 ) 


If Sj^2 is a zero matrix, then and S2J . j ^3 are also all zero matrices. 

The normal use of Eq. 86 is that the plant and model are given and it is 

desired to find z* and u*. 

— — c 

An observer is a form of model following. In an observer, it is desired 

to choose $ stable, F , K , and F so that y = y without using future 
c c c f u ^ 

values of u and using only y , in u* . With these restrictions, Eq. 87 

— ODS — C 

changes to 


K 

F 

c 


0 0 . . ." 


"Sii$ S^^F 0 . .' 

I 

0 

^f. 


c S22 0 . . . 


H D 0 . . 


( 88 ) 


If the equations in Eq. 88 are individually expanded, the standard observer 
conditions are obtained. 


$ S,, 
c 11 


+ 


F C = 
c 


S $ 
11 


( 89 ) 


F 

u 


F S„„ 
c 22 


= ^11^ 


( 90 ) 


40 



K S,, + K.C = H 
c 1 1 f 


( 91 ) 


KfS22 = D 


(92) 


If is applied to Eq. 82, the plant equations become 

•^k+1 ^cllc -^bSjk ^ 

2c.k ' 

and ^observes 


(93) 

(94) 


Comparing Eq. 93 with Eq. 38, the answer to the first question is that 
the compensator is an observer if the compensator matrices satisfy Eqs. 89 
to 92. The answer to the second question is that the compensator/observer 
state, z, , observes S,,x, . 

Equations 89 to 92 imply that if the compensator is to be an observer, 
then the gains in Eq. 93 cannot all be treated as free parameters. For example, 
if and in Eq. 38 are specified, then in Eq. 89 can be computed and 


K is forced to become 
u 

■'u-aT =11'' 


(95) 


Equality constraints, such as Eq. 95, could be used to change the cost 
function in Eq. 48. ‘ 

The answer to the third question can be determined if one compensator is 
treated as the plant and the other compensator is treated as a model so that 
the feedforward matrix equation can be solved. If equality holds in Eq. 87, 

S ^2 ^ 22 ^ are zero matrices and S 22 is an identity matrix, then the compensators 

are alternate internal representations of the same input-output relationship. 

The compensator states are then related by the linear transformation 
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LIST OF SYMBOLS 


Iri general, matrices are represented by capital letters and vectors 
are underscored. 


VARIABLE 


^ii 


w 


a. 

1 


B 

B 


ui 


w 


wi 


b 


1 


EQUATION 

23 

23 

23 

19 

6 

12 

18 

23 

23 

23 

6 

19 

6 

19 

18 


DESCRIPTION 

Reduced-order design wing model system matrix 

Fast actuator model pole 

Reduced-order disturbance pole 

Partitioned, transformed wing model system 
matrix 

52 X 52 wing model system matrix 

Prefilter pole 

Real part of an eigenvalue 

Reduced-order design wing model control matrix 

Control matrix for fast actuator model 

Control matrix for reduced-order disturbance 
model 

52 X 1 wing model control matrix 

Partitioned, transformed wing model control matrix 

52 X 1 wing model disturbance noise matrix 

Partitioned, transformed wing model disturbance 
noise matrix 

Imaginary part of an eigenvalue 
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VARIABLE 

EQUATION 

DESCRIPTION 

C 

2, 28 

Observation matrix 

C 

P 

24 

Reduced-order plant state observation matrix 

C 

w 

7 

1 X 52 wing model state observation matrix 

C . 

Wl 

20 

Partitioned, transformed wing model state 
observation matrix 

D 

22 

_ _-l _ 

D - C „ A„. B „ 
w w2 22 u2 

D 

w 

7 

Wing model control observation matrix 

d(K^) 

58 

Difference between gains 

E 

23 

Disturbance observation matrix 

E 

22 

\ " ^w2 ^22 ®w2 

E 

W 

7 

External disturbance observation matrix 

F 

55 

Positive definite symmetric matrix 

G 

56 

Positive definite symmetric matrix 


76 

Gain schedule gain 


26 

Cross-weighting observation matrix 

^1 

8 

Vertical deflection at node i 

I 

3 

Identify matrix 

J, Jt 

30, 42 


''s’ J 3 

43, 77 

Cost Function 

j 

17b 

cr 

K 

1 

Gain matrix 


36, 37 

Gain matrix from compensator state to plant 



control rate 

Kf, 

31, 36 

Gain from prefilter observation to plant control 
rate 

K 

gs 

76 

Gain scheduled gain matrix 

Sjew 

57 

Projected direction of optimal gain in algorithm 

K 

1 1 

31 

Gain matrix from plant control to compensator 

U1 


control 



VARIABLE 

EQUATION 

DESCRIPTION 

;; K 

V 

31 

Gain from plant control to plant control rate 

K 

y 

31 

Gain matrix from prefilter observation to compen- 
sator control 

$ 

31 

Gain matrix from compensator state to compensator 
control 

k 

1 

Index integer 

L 

9 

Scale of turbulence 

5, 

s 

- 

Number of measurements 

-A, 

M 

29 

Discrete cost function state, control cross- 
weighting matrix 

N 

29 

Number of time steps in cpst function 


77 

Number of design flight conditions 

n 

- 

Number of statps 

P 

49 

Solution mat'pix to cost equation 

/\ 

p 

55 

Intermediate algorithm computation 

Q 

5 

Cross-weighting matrix between plant and 
compensator 

Q 

29 

Discrete cost function state weighting matrix 

Qc 

25 

Weighting matrix for compensator states 


77 

Weighting matrix for gain schedule 


25 

Weighting matrix for wing model states 

q 

- 

Dynamic pressure 

R 

25 

Weighting matrix for wing model control 

A 

R 

29 

Discrete cost function control weighting matrix 

R 

c 

25 

Weighting matrix for compensator controls 

R 

V 

25 

Weighting matrix for wing model control rate 

s 

9 , 

Laplacp transform variable 


86 

Feedforward matrix 
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VARIABLE 


S 

T 

t 

u 

u 

c 

V 

V 

W 

w 

S 


X 




X 




X 

-p 

X 

-qa 

X 

— w 


^ 1 ’ ^2 


EQUATION 

51 

56 

17a 


1 

3 

45b 

9 

13 

44b 

80 

27 

47 


27 

8 

8 

8 

23 

80 

23 

8 

6 

19, 20 


DESCRIPTION 

Solution matrix to covariance like 
equation 

Intermediate algorithm computation 

Block diagonalizing transformation matrix 

Time 

Control State 
Compensator control state 
Measurement noise covariance 
Aircraft velocity 
Control rate 

Plant process noise covariance 

Diagonal weighting matrix in gain schedule 

Discrete plant prpcess noise 

Initial condition covariance 

Control actuator state initial condition 
covariance matrix 

State vector 

Actuator state 

Disturbance state vector 

Dryden wind model state vector 

Reduced and balanced disturbance/gust state 
vector 

Complex eigenvector 

Plant state vector for reduced model 

Unsteady lift and moment state vector 

52 X 1 wing model state vector 

Partitioned, transformed wing model state 
vectors 
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VARIABLE 

EQUATION 

DESCRIPTION 

— 

yobs*y 

85b, 2 

Measurement vector 



4 

Compensator state measurement vector 


yf 

12 

Prefilter state measurement 

— 


3 

Compensator state vector 


VARIABLE 

(GREEK) 

EQUATION 

DESCRIPTION 



a. 

1 

59 

Algorithm stability parameter 



3 

13 

Continuous model compensator stability 
parameter 

' 

^1’ ^2 

3 

Discrete model compensator stability 
parameter 

I 

r 

27 

Plant discrete control matrix 

— 

r 

c 

82 

Compensator control matrix 


r 

p 

27 

Wing model discrete control matrix 


r 

u 

35 

Control discrete control matrix 

, 

r 

w 

27 

Wing model discrete disturbance matrix 


Y 

2 

Measurement noise 

— 


7 

Accelerometer Measurement Noise 



8 

Actuator position 


\ 

6 

Actuator housing white noise disturbance 
state 

p 

X 

16 

Complex eigenvalue 


00 

9 

Gust state vector 


I 

29 

Summation 


a 

- 

Number of compensator states 


a 

Wg 

9 

Variance of gust 


e 

53 

Small number greater than zero 
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VARIABLE 

(GREEK) 


$ 

$ 


c 



$ 


P 

$ 


u 


At 


EQUATION 

27 

82 

50 

27 

33 


DESCRIPTION 

Plant state transition matrix 
Compensator state transition matrix 
Closed-loop plant matrix 
Wing model state transition matrix 
Control state transition matrix 
Control and measurement sample interval 


SPECIAL - OPERATORS 


DESCRIPTION 


(•) 

(_) 

8( )/a( ) 

( )* 

( )* 

00 

I 

e[ ] 


Derivative of quantity with respect to time 
Vector quantity 

Partial derivative of one variable with respect 
to another 

Complex conjugate transpose 

Star trajectory 

Infinity 

Integral 

Expected Value 

Transpose 

Inverse 


ACRONYM 

LQG 

LQR 

tr 


CORRESSPONDING PHRASE 
Linear quadratic Gaussian 
Linear quadratic regulator 
Trace of a square matrix 
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TABLE 1 POLES AND ZEROES OF THE WING MODEL 



OPEN-LOOP EVALUATION 
MODEL POLES 

EVALUATION MODEL 
ZEROES 





10.8 ± j59.8 

13.9 


-40.0 ± j78.5 

-50.0 


MODES USED ■ -5.8 ± J162.3 

-4.2 + J163.0 * 


FOR BODE • -40.0 ± j225.5 

-39.4 + J226.9 * 


plot -125.0 ± j352.3 

-228.0 ± J1019.0 


-26.4 ± j377.3 

-21.9 ± J370.0 


-16.5 ± j418.2 

-5.1 ± J415.2 


-37.4 ± J557.0 

-19.0 ± J545.0 






-36.6 ± j663.0 

-36.3 ± J659.0 * 


-34.6 ± J776.0 

-4.4 + J770.0 


MODES THAT -39.7 ± j 820.0 

-54.7 ± J812.0 


FOLD FOR -59.0 ± j 908.0 

-59.0 ± j906.0 * 


At = 0.005 -16.0 ± J1278.0 

-89.7 ± jl338,0 

MODES 

-17.7 + J1862.0 

-32.6 ± J1891.0 

RESIDUALIZED 

-20.8 ± J2638.0 

-21.2 ± J2639.0 * 

AND not USED 

-28.4 + J5617.0 

0.0 ± jO.O 

IN REDUCED 



MODEL DESIGN 

. 


-228.0 

-228.0 * 


-232.0 

-232.0 * 


-254.0 

-254.0 * 


-266.0 

-265.0 * 


-284.0 

-284.0 * 


-294.0 

-294.0 * 


-320.0 

-321.0 * 


-332.0 

-332.0 * 


-364.0 

-369.0 * 


-379.0 

-381.0 * 


-421.0 

-482.0 


-426.5 

-425.0 * 


-501.7 ± j2.3 

-497.0 


-706.0 

-714.0 ± J34.9 


* The pole and zero almost cancel for the actuator 
to accelerometer output wing model transfer 
function. 
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TABLE 3 OPEN AND MAPPED CLOSED-LOOP POLE LOCATIONS 


. 

OPEN^LOOP EVALUATION 
MODEL POLES 

MAPPED CLOSED-LOOP EVAL- 
UATION WING MODEL . POLES 

MAPPED CLOSED- 
LOOP REDUCED-ORDER 
WING MODEjL POLES 

'^ODE^SED 10.8 ± j59.8 

-14.8 ± J59.7 

-15.5 ± J59.7 

IN REDUCED -40.0 ± J78.5 

-76.8 ± jl39.0 

-79.2 ± jl43.0 

MODEL DESIGN -5.8 ± j 162.3 

-6.8 ± J161.0 

-6.9 ± jl61.0 

-40.0 ± J225.5 , 

-41.0 ± j225.0 

-41.0 ± J225.0 

^ -125.0 ± J352.3 

-149.0 ± j442.0 

-136,0 ± J435.0 

MODES USED -26.4 ± j337.3 

-26.0 ± J375.0 


FOR BODE -16.5 ± J418.2 

-15,4 ± J416.0 


PLOT -37.4 ± j557.0 

-38.8 ± j558.0 





MODES THAT”^ r36.6 + j 663,0 

-36.4 ± J593.0 


FOLD FOR -34.6 ± J776.0 

-34.6 ± j480.0 ' 


At = 0.005 -39.7 ± J820.0 

t39.7 ± J436.0 


-59.0 ± j908.0 

-59.0 ± J348.0 


-16.0 ± J1278.0 • 

-15.7 + j21,8 


-17.7 ± jl862.'0 

-17.7 ± j606.0 


-20.8 ± J2638.0 

-20.8 ± J125.0 


-28.4 ± J5617.0 

-28.4 ± J591.0 





-228.0 

-228.0 


-232.0 

-232.0 


-254.0 

- 254.0 


-266.0 

-266.0 


-284.0 

-284.0 


-294.0 

-294.0 


-320.0 

-321.0 


-332.0 

-332.0 


-364.0 

-348.0 


-379.0 

-381.0 


-421.0 

-538.0 


-426.5 

-425.0 


-501.7 + j2.3 

-t497.0, -485.0 


-706.0 

-725.0 
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TABLE 4 NQNZERO DIAGONAL ELEMENT DESIGN PARAMETERS 




WEipHT 

COVARIANCE 

Wing Mpde 1 

1 

I-* 

(2.24)^ 


Wing Mode 1 

-^o 


(3.16)^ 

4 Coippensator 
States 

- Qc 

(0.316)^1 


4 Compensator 
States 

- X 

o 


(10. 0^)1 

Actuator Control 
Position 

- R 

(1.0)^ 


Compensator 

Control Position 

- R 

c 

(l.O)^I 


Actuator Control 
Rate 

- R 

V 

(1.23)^ 


Actuator Control 
Position 

Accelerometer 

Measurement Noise 

- X 

u 

- V 

a 


TABLES 5,7, 
8,9,10,11,12 

(0.1)^ 

Cross Coupling 
Weigl^t 

- Q 

(31.6)^ 





■■“1 - '1 '1 1 ! ^ 1 ] ■ ■] 1 1 '1 ] 1 1 1 1 ■ 1 ■ 1 


TABLE 5 EFFECT OF TARYING INPUT NOISE COVARIANCE 


INPUT NOISE 
COVARIANCE 

■ EVALUATION MODEL 
RMS RESPONSE 

GAIN MARGIN 
REDUCED MODEL 

PHASE MARGIN 
REDUCED MODEL 

BANDWIDTH . 
0 - 6dB 

Xu 

unns 
deg ’ 


*^rins* 

deg 

• 

deg7sfec 

dB 

rad /sec 

dB 

rad/ sec 

deg 

rad/sec 

deg 

rad/sec 

rad/ sec 

0.0 

7.0 

380 

4.2 

233 

2.4 

61 

-4.2 

98 

-13.3 

51 

28.4 

78 

107 

0.1 

10.0 

392 

5.8 

227 

5.0 

60 

-7.5 

156 

-27.5 

46 

50.6 

81 

135 

1.0 

11.4 

392 

6.6 

227 

5.6 

60 

-8.0 

190 

-34.0 

45 

55 . 0 

81 

154 


\-n 

Ln 




TABLE 6 VARIANCE OF NODE DEFLECTIONS WITH 0.3048 m/s 
WIND GUST FOR THE NOMINAL CONTROL DESIGN 


1 

NODE 

2 

POSITION (CENTIMETERS) 
3 4 5 6 

7 

0.013 

0.10 

0.25 0.47 

0.74 1.1 

1.4 

1 

NODE ROTATION 
2 3 4 

(DEG) 

5 6 

7 

0.02 

0.08 

0.13 0.19 

0.27 0.38 

0.48 




1 


1 


1 


1 


■] 


1 


TABLE 7 EFFECT OF LOWERING SAMPLE RATE 


STARTING 

GAIN 

DIFFERENT 

SAMPLES 
PER SECOND 

INPUT 

NOISE 

COVARIANCE 

X 

u 

EVALUATION MODEL 
RMS RESPONSE 

GAIN MARGIN, 
REDUCED MODEL 

PHASE MARGIN, 
REDUCED MODEL 

BANDWIDTH 
@ - 6dB 

1.0/At 

deg 

deg/ sec 

dB 

rad/ 

sec 

dB 

rad/ 

sec 

deg 

rad/ 

sec 

deg 

rad/ 

sec 

rad/sec 

NO 

200 

0.1 

5.8 

227 

5.0 

60 

-7.5 

156 

-27.5 

45 

50.6 

81 

142 

NO 

160 

0.25 

5.9 

219 

5.1 

60 

-6.3 

155 

-30.5 

46 

48.5 

80 

151 

YES 

140 

0.25 

5.5 

210 

5.0 

60 

-6.0 

152 

-30.2 

47 

45.0 

80 

148 

YES 

100 

0.25 

5.1 

171 

4.4 

60 

-5.5 

138 

-29.0 

48 

37.0 

78 

NONE 


Ln 




TABLE 8 EFFECT OF PREFILTER POLE VARIATION 


STARTING 

GAIN 

DIFFERENT 

PREFILTER 

POLE, 

INPUT 

NOISE 

COVARIANCE 

Xu 

EVALUATION MODEL 
RMS RESPONSE 

GAIN MARGIN 
REDUCED MODEL 

PHASE MARGIN, 
REDUCED MODEL 

BANDWIDTH 
@ - 6dB 

-1 

sec 

'^ms’ 

deg 

6 

des7sec 

dB 

rad/ 

sec 

dB 

rad/ 

sec 

• 

deg 

rad/ 

sec 

deg 

rad/ 

sec 

rad/sec 

NO 

-60.0 

0.1 

5.4 

222 

4.7 

59 

-8.2 

149 

-26.0 

46 

45.6 

79 

120 

NO 

-50.0 

0.1 

5.8 

227 

5.0 

60 

-7.5 

156 

-27.5 

45 

50.6 

81 

X42 

NO 

-40.0 

0.1 

5.9 

226 

5.1 

58 

-7.3 

156 

-28.0 

45 

51.9 

81 

143 

YES 

-30.0 

0.1 

6.1 

- 

225 

5.4 

58 

-7.1 

157 

-30.0 

46 

52.1 

82 

144 




1 1 ■ 1 1 1 1 'I 1 1 1 1 1 ] "1 1 1 1 1 “] 


TABLE 9 EFFECT OF DIFFERENT CONTROL STRUCTURE 


STARTING 

GAIN 

DIFFERENT 

CONTROL 

STRUCTURE 

INPUT 

NOISE 

COVARIANCE 

EVALUATION MODEL 
RMS RESPONSE 

GAIN MARGIN, 
REDUCED MODEL 

PHASE MARGIN, 
REDUCED MODEL 

BANDWIDTH 
@ - 6dB 

^rms’ 

deg 

^rms> 
deg /sec 

dB 

rad/ 

sec 

dB 

rad/ 

sec 

deg 


deg 

rad/ 

sec 

rad/sec 



\ 

NO 

FIG. 2 

0.1 

5.8 

227 

5.0 


60 

1 

156 

-27.5 

45 

50.6 

81 

142 

NO 

FIG. 6 

0.01 

5.6 

228 

5.0 


57 

-7.8 

136 

-24.5 

46 

44 

87 

120 


Ln 

VO 



























o> 

o 


TABLE 10 CONSTANT GAIN DESIGN WITH VARYING DYNAMIC PRESSURE 


DYNAMIC 
PRESSURE 
q, kPa 

INPUT 

NOISE 

COVARIANCE 

EVALUATION MODEL 
RMS RESPONSE 

GAIN MARGIN, 
REDUCED MODEL 

PHASE MARGIN, 
REDUCED MODEL 

BANDWIDTH 
@ - 6 dB 

^rms , 
deg 

^rms , 
deg/sec 

dB 

rad/ 

sec 

dB 

rad/ 

sec 

deg 

rad/ 

sec 

deg 

rad/ 
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rad/sec 

X 

u 
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TABLE 11 OPTIMAL OUTPUT FEEDBACK DESIGNS FOR VARYING DYNAMIC PRESSURES 


STARTING 

GAIN 

DIFFERENT 

DYNAMIC 
PRESSURE 
q, kPa 

INPUT 

NOISE 

COVARIANCE 

EVALUATION MODEL 
RMS RESPONSE 

i;ain margin, 

REDUCED MODEL 

PHASE MARGIN, 
REDUCED MODEL 

BANDWIDTH 
@ - 6 dB 

^rms , 
deg 

*5rms , 
deg/sec 

dB 

rad/ 

sec 

dB 

rad/ 

sec 

deg 

rad/ 

sec 

deg 

rad/ 

sec 

rad /sec 

X 

u 

NO 

8 
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TABLE 12 SCHEDULED GAIN DESIGN WITH VARYING DYNAMIC PRESSURE 


DYNAMIC 
PRESSURE 
q, kPa 
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EVALUATION MODEL 
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^rms, 

deg 
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FIGURE 2 CONTROL LAW BLOCK DIAGRAM 
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FIGURE 8 OPEN-LOOP DYNAMIC PP^SSURE ROOT LOCUS AT MACH 0.9 
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FIGURE 9 CLOSED-LOOP DYNAMIC PRESSURE ROOT LOCUS FOR CONSTANT GAIN DESIGN 
AT q = 8,0 kPas 
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